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Abstract 

A theory of mesoscopic fluctuations in disordered thin superconducting films 
in a parallel magnetic field is developed. At zero temperature and at a suf- 
ficiently strong magnetic field, the superconducting state undergoes a phase 
transition into a state characterized by superfluid densities of random signs, 
instead of a spin polarized disordered Fermi liquid phase. Consequently, in 
this regime, random supercurrents are spontaneously created in the ground 
state of the system, which belongs to the same universality class as the two 
dimensional XY spin glass. As the magnetic field increases further, meso- 
scopic pairing states are nucleated in an otherwise homogeneous spin polar- 
ized disordered Fermi liquid. The statistics of these pairing states is universal 

depending on the sheet conductance of the 2D film. 
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I. INTRODUCTION 

Recent experiments on thin superconducting films in a parallel magnetic field |l|] have 
rekindled interest in this field. If the thickness of the films d is small enough, the orbital 
effect of the magnetic field can be neglected and the suppression of superconductivity in the 
film is due to the Zeeman effect [0,^3 • It has been observed that the resistance of such films 
at low temperatures and high enough magnetic fields exhibits very slow relaxation in time 
||1||. This behavior is characteristic for spin and superconducting glasses. Below we discuss 
a possibility that mesoscopic fiuctuations of superconducting parameters in disordered films 
account for such a behavior. 

Mesoscopic physics in a noninteracting electron system has been known for a while P,§ . 
The energy spectrum in a mesoscopic sample was shown to exhibit Wigner-Dyson statistics, 
which is universal, only dependent on the symmetry of the Hamiltonian 0. The long range 
level repulsion in the energy spectrum leads to a suppression of fiuctuations of levels within 
an energy band of width Ec = -D/-L^(Thouless energy). L is the length of the sample, 
D = vpl/3 is the diffusion constant of the film, vp is the Fermi velocity, I is the elastic mean 
free path. For an open sample, the fiuctuation of number of levels 6N within the energy 
band of Thouless energy Ec is of order unity, 

for a 2D film, where the corresponding average number of levels < N >= E^dvoEc = G, 
with z/q being the average density of states in the metal on the Fermi surface. /3 is a factor of 
order unity depending on the symmetry of the Hamiltonian. G = kpdl, is the dimensionless 
conductance of the 2D normal metal in units of e^/h, kp is Fermi wave length and the brack- 
ets denote averaging over realizations of random potential. Consequently, the transport 
is governed by UCF (universal conductance fiuctuation) theory. The conductance exhibits 
sample specific fluctuations, with amplitude e^/h, independent of the average conductance 
of the sample |^J§]. More generally, any physical quantity in a mesoscopic sample consist of 
an ensemble average part and a sample specific part due to quantum interference. 



On the other hand, disordered superconductors have been studied long ago 0. It was 
shown that the ground state condensate wave function is homogeneous and the critical 
temperature remains unchanged in the presence of weak nonmagnetic disorders. To derive 
the dirty superconductor theory, one has to assume that 1). the effective interaction constant 
in the Cooperon channel remains the same as in a clean superconductor; 2). the condensate 
wave function is translationally invariant; 3). the time reversal symmetry is preserved. The 
first assumption, though is not true in the thin film limit where the Coulomb interaction in 
Cooperon channel can be greatly enhanced, is valid in the bulk limit [§ . We will assume its 
validity because it does not affect the result present in this paper as far as the renormalized 
interaction constant is still negative. The translation invariance is not a generic symmetry of 
the original Hamiltonian in the presence of disorder and the second assumption is true only 
after the impurity average is taken in the semiclassical limit. The sample specific quantum 
interference effect which is of the same origin of Wigner-Dyson statistics was not taken into 
account. The consequency of such an effect which breaks the translation invariance is one 
of the subjects of this article. 

The most fundamental aspect of Anderson theory for a dirty superconductor is the 
absence of spontaneous time reversal symmetry breaking] that is, the stability of a BCS state 
with respect to possible frustrations, even when / is much shorter than the coherence length. 
This is in contrast to how a ferromagnet responses to impurities: A spin glass phase which 
does not have a conventional long range order but does have Edward-Anderson type long 
range order always takes over when impurities are added into the system [^]. This central 
issue will be addressed in this article, in connection with nodes in the spatial dependence of 
exchange interactions and the distribution function of the exchange interactions. 

Unlike in the noninteracting metal where the mesoscopic physics is relevant only in a 
finite sample smaller than the dephasing length, in the presence of off-diagonal long range 
order, it reveals itself in the thermodynamic limit. However, when the elastic mean free path 
/ exceeds the Fermi wave length h/kp, mesoscopic fluctuations of various physical parameters 
of superconductors are smaller than their averages []l0|-[r3[| . Thus, it seems that they hardly 



affect macroscopic observable quantities. It was realized later that there are situations where 
mesoscopic ffuctuations determine macroscopic properties of a superconducting sample. One 
example is a superconductor in a magnetic field close to the upper critical field ifc2, where 
the magnetic field dependence of the superconducting critical temperature is determined by 
the mesoscopic fluctuations |T^. In general, the mesoscopic effects are not only relevant in 



a disordered superconductor but also determinant to the global phase rigidity. 

In this paper we consider the case, where the magnetic field is parallel to the thin 
superconducting film and the main contribution to the suppression of superconductivity by 
the magnetic field is due to Zeeman splitting of electron spin energy levels. We show that 
at low temperatures T and high enough magnetic fields H, parallel to the film, the system 
exhibits a transition into a state where the local superfluid density Ns{v) (which is the ratio 
between the supercurrent density J^ and the superfluid velocity V^) has a random sign. In 
this case the system belongs to the same universality class as the two-dimensional XY spin 
glass model with exchange interaction of random signs. We also find that as the magnetic 
field is decreased from above the critical field, mesoscopic pairing states are nucleated in an 
otherwise spin polarized disordered Fermi liquid. The characteristic length scale at which 
pairing takes place increases as the critical field is approached. The statistics of these pairing 
state is universal depending on the sheet conductance only. 

The idea that the superfluid density can be of random signs has a long history 
T5| , p!0| - |T2| , p!6| -[T8[| . However, in the absence of magnetic fields and at zero temperature in 



disordered superconductors (^o ^ ^ ^ Ti/kp) the variance of the superfluid density, aver- 



aged over the superconducting coherence length ^o = JD/Aq, turns out to be much smaller 
than its average [p!0HT3|| 

^<^I^-±«l (2) 

Aq is the value of the order parameter ai T,H = 0. Similarity between Eqs.1,2 suggests 
the intimate relationship between the fluctuations of the superfluid densities and universal 
Wigner-Dyson statistics. In fact, Eq.2. follows as a consequency of the fluctuation of number 



of levels within Thouless energy band in a volume of size of the coherence length ^o- As 
long as kpl ^ h, the regions where the superfluid density is negative are rare and do not 
contribute significantly to macroscopic properties of superconductors. The situation in the 
presence of a magnetic field parallel to the film is different, because the average superfluid 
density decays with H faster than its variance. Hence, at high enough magnetic field the 
amplitude of the mesoscopic fluctuations of Ns{r) becomes larger than the average, and the 
respective probabilities of having positive and negative signs of Ns{r) are of the same order 
even at kpl/Ti ^ l(See below). This was first pointed out in an early paper by the author 

In section 2, we present the qualitative picture of this phenomenon, emphasising on the 
sensitivity of mesoscopic fluctuations of spin polarization energy to the change of the pair 
potential. In section 3, we study the mesoscopic fluctuations of the superconducting order 
parameter near the critical regime and show that there are spontaneously created currents 
in the ground state. In section 4, we derive the distribution function of the ground state 
condensate wave function at a magnetic field higher than the critical one. In section 5, 
we discuss the role of the exchange interaction. In section 6, we discuss the mesoscopic 
effects in a finite size superconductor. In conclusion, we propose possible experiments to 
observe these effects and point out a few open questions, including implications on d-wave 
superconductors. 

II. QUALITATIVE PICTURE 

A theory of magnetic field induced phase transition which does not take into account 
mesoscopic fluctuations predicts pO| , pT| ,^ that at low temperatures the superconductor- 



normal metal transition, is of first or second order depending on whether the parameter 
AqTso is larger or smaller than unity respectively. Here, Tso is the spin-orbit relaxation 
time. The if-dependence of the order parameter for the two limiting cases was discussed 
extensively in p]| , pT| . 



From now on, we restrict ourselves to the limit Aot^o ^ 1, where the theory predicts a 
second order phase transition between the superconducting state and the normal state. At 
T = and within an approximation which neglects mesoscopic effects, the value of the crit- 
ical magnetic field if ° is the result of the competition between the average superconducting 
condensation energy density (-Ec) ~ J^o^o ^^^ ^^^ polarization energy of the electron gas 
in the magnetic field. The average spin polarization energy density of nonsuperconducting 
electron gas is of order {Ep{0)) ~ uq^hbHY- Its relative change in the superconducting 
state is of order (^p(O)) - {Ep{A)) ~ |7rAor,o(Ep(0)) < (f^p(O)) ^^. As a resuh we get 



an expression for the critical magnetic field H^ = Hcc^^qTso)^"^ ^ Hcc- Here Hcc = Ao/fiB 
is the Chandrasekar-Clogston critical magnetic field of the superconductor-normal metal 
transition for AqT^q -^ oo and fiB is the Bohr magneton. 

Consider the mesoscopic fluctuations of the quantities, discussed above, in a volume 
whose size is of the order of the coherence length ^o- To calculate the amplitude of mesoscopic 
fluctuations of the polarization energy 6Ep, we use the conventional diagram technique for 
averaging over realizations of random potentials [^. By evaluating the diagrams in Fig.2.d 
(see Appendix A), we have 

\/((^Woc^(Aor,„)(E,(0)). (3) 

This part of the polarization energy is sensitive to the change of the pair potential (if ~ He) 
just as the quantum interference effect is sensitive to the change of impurity potentials. The 
A-dependent part of the mesoscopic fluctuation of spin polarization energy can be obtained 
by calculating the diagrams in Fig.2.d, 

1 , . . 5A' 



{{6Ep{6A) - 6Ep{A = O))^) ex -(r,,Ao) < E^iO) > —^. (4) 

Linear in 6 A term vanishes because a quasiparticle is reflected into a quasi hole when it is 
scattered by 6 A. As a result, the change of mesoscopic fluctuations of the spin polarization 
energy associated with the change of the pair potential 6 A of order A{H) is given as 

{{6Ep{A{H)) - 6Ep{A = O))^) ex ^(('^^^(O))^) (^) • (5) 



Eq.5 shows that below the critical field, though the average cost of condensation energy 
and kinetic energy to have a configuration with order parameter equal to zero in some 
regions are positive, the mesoscopic fiuctuations of polarization energy associated with such 
a configuration are of random signs. Since both the polarization energy and the condensation 
energy are fiuctuating quantities, A(r) should also be spatially fiuctuating. Particularly, 
when H^ is approached, the average cost of energy vanishes, the spatial structure of the 
order parameter is completely determined by mesoscopic fiuctuations of spin polarization 
energy. One elaborates the following argument to confirm this picture. 

Consider a domain of size Ld ^ ,^0 where the value of A(r) differs from its bulk value 
by a factor of order unity. An estimate for the energy of such a domain consists of three 
terms, namely 

SE{A) = uoAiHfdiC^Ln^o + C^^^^^L^ + ^3^0') (6) 

where Ci, C2, C3 are factors of order of unity. The first term in Eq.6 corresponds to the A- 
dependence of mesoscopic fiuctuations of polarization energy and has a random sign. When 
estimating this term we have taken into account that regions of size ^0 make independent 
random contributions into Eq.6. The second and third term are the average condensation 
energy and surface (gradient) energy of the domain, respectively. It follows from Eq.6 that 
when Lo ~ ^{H) = C,oA.o/{A{H)), there is an interval of magnetic fields near the critical 
one 

Hc-Hr^^ (7) 

where the first term is larger than the second and the third ones. Here (A (if)) = 



A()J(H^ — H)/H^ is the average superconducting order parameter. It means that, in this 
case the spatial distribution of A(r) is highly inhomogeneous and the amplitude of the 
spatial fiuctuations of A(r) is of order of its average, while the characteristic size of the 
domains is of order of Ld ~ ^{H = H^{^ — ^/G^)) ~ ^oG. At H far away from H^, the 
typical mesoscopic fiuctuations are smaller than the average contribution and do not change 



the most probable configuration; they only introduce exponentially small concentrations of 
defects originating from the statistically rare events. 

The second kind of the instability which is intimately connected with this inhomogeneous 
mesoscopic superconducting state is spontaneous creation of long range current. In other 
words, the superfluid density, which is the second derivative of the energy with respect to 
superfiuid velocity in this region has a random sign and is no longer positive defined. To see 
this, one should consider states with finite superfiuid velocity V^ = ^(Vx + 2e/cA), where 
x(r) is the phase of the order parameter, m is the electron mass and A(r) is the vector 
potential, which has a direction parallel to the film and m is the electron mass. If Vs(r) is 
of the order of the critical velocity, all three terms in Eq.6 are modified by factors of order 
of unity when compared with the case V^ = 0. The second and the third term in Eq.6 
decrease with V^, while the first term is changed in a random direction. This means that at 
high enough magnetic fields, states with nonvanishing value of Vs(r) can have lower energy 
than the states with V^ = 0, and that the system is unstable with respect to the creation 
of supercurrents flowing in random directions. In this estimate we neglect the energy of 
the magnetic field associated with Vs(r) in the thin film limit. Since at each point of the 
system the possible energy gain associated with finite value of Vs(r) is independent of the 
direction of V^, the ground state of the system is highly degenerate and belongs to the same 
universality class as XY spin glass with a random sign of the exchange interaction. 

In principle, both the spin polarization energy and the condensate energy fluctuate from 
region to region. Both are related to the density of states, though the condensate energy 
should be determined by the solution of the self consistent equation. For the above argument 
to be true, we have to assume that the spin polarization energy and the condensate energy 
fluctuate independently. The argument can be carried out in a similar fashion even if these 
two are partially correlated, as far as they are not fully correlated [^ . More serious consid- 
eration of the existence of inhomogeneous state is addressed in term of the self consistent 
equation in the next section. 

It is important to mention that even in the case of small magnetic flelds in the presence of 



spin orbit scattering the time reversal symmetry is broken and the electron wave functions 
are complex. Therefore there are currents in the ground state of the system which have 
random directions. These currents exist even in normal metals. Diagrammatic calculation 
leads to an expression for the correlation function of the current density J(r) in a normal 
metal induced by a magnetic field (|r — r'| ^ /i/Zcj?) as shown in [|r^ 

(J,(r)J,(r')) ^ 5,,^^tt,,{^^bH)%v - r'). (8) 

Here r = l/vp is the elastic mean free time. It is important to note, however, that for a 
given configuration of the scattering potential and at a given value of the external field the 
spatial distribution of J(r) is a unique function. This implies that the currents described 
by Eq. 8 do not exhibit features which can be associated with superconducting glass states. 
In other words, at small H the superfiuid density of the superconducting state is positive 
which means that states with nonzero superfiuid velocity have larger energy than the ground 
state. The rare regions, where A^s(r) < and the supercurrents exist in the ground state, 
are screened effectively due to the Meissner effect. They do not affect significantly the 
macroscopic behavior of a sample. 

III. CONDENSATE WAVE FUNCTION: H < H^ 

Below, we will be interested in supercurrents much larger than those described by Eq.8. 
Such currents are spontaneously created at strong enough magnetic fields as a result of the 
instability associated with the random sign of superfiuid density. To evaluate the variance 



of the superfiuid density we consider the Gorkov equation for A(r) p5|] , 

A{r) = g I dr'Kir,r';H,Air),{A{r)})A{v'), (9) 

where g is the dimensionless interaction constant, 

K{r,r';H,A{r),{A{r)}) 

= u^'kT Y: G^'^r, r'; H, A(r), {A(r)})a,"^"3G'^3«4(r, r'; H, A(r), O^^^S (10) 



G°^(r, r'; iJ, A(r), A(r)) is the exact one particle electron Matsubara Green function in 
the presence of pair potential A(r), a and /5 are spin indexes, a^^ is the y component of 
Pauli matrix and e = (2n + V)TxkT is the Matsubara frequency. Both A(r) and Kiv^ r') in 
Eqs.9 and 10 are random functions of the realizations of scattering potential in the sample. 
Averaging Eq.lO over realizations of the random potential and using the approximation 
(A(r)A'(r, r')) = {A{r; H)){K{r,r',H)) we get the above mentioned expression for H^. 

In the case of strong magnetic fields, when A{H, r) ^ Aq, we can expand Eq.lO in terms 
of A(r). Since A(r) varies slowly over distances of the order of ^O; while {K{r,r')) decays 
exponentially for |r — r'| ^ ,^0, we can also make the gradient expansion of Eq.lO. As a 
result we get from Eq.lO 

(eo^(V - .^A)^ + ^^^A(r)) + / 6K%r, r', H, A) A(rOrfr' = ^, (11) 

where 

<5ir°(r, r') = ir°(r, r') - (ir°(r, r')) (12) 

and K^{r,r') = K{A{r) = 0,r,r'). The difference between Eq.ll and the conventional 
Ginsburg-Landau equation is the third term in Eq.ll which accounts for mesoscopic fluctu- 
ations of the kernel i^°(r, r'). It is precisely this term, which at high magnetic flelds leads 
to the random sign of superfluid density. To proceed further, we calculate the correlation 
function 

C(ri, r[; r^, r^) = {6K\r,, r[)6K\r,, r^)) (13) 

using the diagrams shown in Fig.2f. And we obtain 

C{ri,r[;r2,r'^) = (27rT)^— ^exp(iq ■ ri - iq_i -rs + iq'i ■ r'^ -ic[r[) 



5(q + q'l - q' - qi)5(Qi - Q2 + q " q') 



a\-a\ /3i -/3i ft -ft _Q!2 -^2 



^^(Z^qi • q + 2|ei + t^^D^ ■ q^ + 2|ei + e2|)<^-"V^^^-^V^^ '-cy^ 

\c^'' (q)c^2^r"^ {'^)CIT (qOC^.?''^ (q'lOc^J-t? (Qi)C7;i? (Q2) 
^c^r^ (q)C2^r^ (q')C2T^ {'^r)cit (q'i')^JI? (Qi)/^j;? (Q2)] . (m) 
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The large distance asymptotics of the correlation function in Eq. 14 takes the form, 



C(ri,r;;r2,r2) 



lecs 



eoS{r^ - r2)S{r[ - r',)6ir^ - r[) + 5(ri - r',)Sir[ - v,) 



|ri-r'J4 



■ (15) 



Eq.l5 characterizes the sample specific interference effect on the Cooperon propagator 
defined in Eq. 10. It determines the mesoscopic fluctuations of the superconducting order 
parameter, which represent the deviation of the exact ground state from the translation- 
ally invariant state. Employing the perturbation theory with respect to 5K^{Y^r') we get 
from Eq.ll an expression for the correlation function of the mesoscopic fluctuations of the 
superconducting order parameter (5A(r; H) = A(r; H) — {A{H)) 

{5A{r)5A{r')) = A\H) j dv^dr'^dv^dv'^Kv, ri)L(r', r;)C(ri, r'^; rs, r^) 

Hrr^-V ^^PN(^-ri)) Hfil 

Taking into account Eq.l5, 



A| J 1 - i Hi{H)/i,)[\v - v'\/i{H)f, |r - r'l « ^H); 



(5A(r)5A(r')) oc ^ ^ 2 v v .. /u ^^^^ 

exp[-|r-r'|/e(i/)],|r-r'|»e(i/) 



It follows from Eq.l7 that the amplitude of the fluctuations of the order parameter in the 
two-dimensional case is almost independent of iJ, but the average order parameter decreases 



with H. As a result, perturbation theory holds as long as {A{H)) / Aq = J{H^ — H)/H^ ^ 
G^^. This also justifies Anderson theory of dirty superconductors in the absence of an 
external magnetic field: the ground state is approximately translationally invariant though 
the translation invariance is not a generic symmetry of original Hamiltonian in the presence 
of impurity potentials. 

Eq.l7 implies that a homogeneous superconducting state becomes unstable against the 
mesoscopic fiuctuations near the critical point. Such an instability against the inhomoge- 
neous state can also be visualized if the magnetic field is decreased from above the critical 
field. The generalized curvature characterizing the stability of a metal or at A(r) = 0, is 
defined as, 

11 



_ ^'g({A(r)}) 
^ • ' iA(r)«A(r') ^ ' 

where E{{A}) is the energy of a configuration {A(r)}, and 6/6A{r) stands a functional 
derivative. Following Eq. 11, we obtain 

0{r, t') = -(eo^(V - t^Af + ^^^^)5(r - r') - 5K'{t, r') (19) 

Eq.l9 shows that 

< {0{r,,T[)- < 0{r,,r[) >){0{r,,r',)- < 0{r,,r',) >) >= C(ri, r',; r^, r'^). (20) 

Generally speaking, the curvature matrix in Eq. 18 is not a positive defined one because of 
fiuctuations. However, well above the critical point, the probability to find the region where 
the curvature is negative is exponentially small; the ground state will be a normal metal 
with exponentially small concentration of superconducting droplets. The probability to find 
these droplets will be discussed in detail in the next section. Here we want to point out that 
when the critical point is approached, the mesoscopic fiuctuations of the curvature matrix 
becomes larger than the positive defined average part (first two terms) and the probability 
of finding the superconducting regions in the ground state becomes of order of one. This 
again implies that the most probable configuration near He is an inhomogeneous state. 

The instability against spontaneous creation of current state can be demonstrated via 
studying the superfiuid density defined as 

which can be written in term of exact Green functions when A(iJ) ^ Aq 

e^kT 



Ct2—OL2 

y 



m ^ J ^ 

^^^^"^(rg, r')Vr'G'2^"^(r', r^)A{v,)a'^^~^'GzZ'''' (r4, r). (22) 

Expanding Eq.22 in terms of 6A{r,H) <^ A{H) we get an expression for the nonlocal 
superfiuid density A/'*''(r, r'), which is valid as long as {A{H)) / Aq ^ G^^ and 6Ns = {Ng — 

(iV,)) « (iV,) 

12 



M'^iv, v') =< Af'Jir, r') > +6Af'J{r, r'). (23) 

The average superfluid density is (5-correlated over distances larger than /, 

(AC^-(r, r')) = N^^^S,,6{r - r'), (24) 

where iV° = eN{l/C,o)'^ is the average superfluid density at if = and A^ is the electron 
concentration in the metal. 

(5A/ij(r,r') is determined by the diagrams in Fig.2g,2h. Following the calculations in 
Appendix B, 

{6^Hr„ r[)SArp'{r2, r',)) (5A(ri)5A(r2))2 



{N,{H)y {A{H)Y 



-6ij6iff6{ri - r[)6{r2 - r'^ 



1 t4 

+ Tfi?^ I ^^Tui^ii'^j'j'^i^i - ^2)5{r'-^ - r'a) + 5ij>5i'j5{ri - T'2)5{r[ - ra)). (25) 

Here Ns{H) = N^ {A{H) / Aq)"^ is the superfluid density at H. The flrst term in Eq.25 
is connected to the fluctuations of the order parameter A(r) in Eq.l6 (the corresponding 
diagram is shown in Fig.2g). The second term in Eq.25 is related to the fluctuations of 
the Green functions Gi^{r,r') (the corresponding diagram is shown in Fig.2h). When the 
magnetic fleld is close to the critical one, i.e. \H — H^\/H^ ~ G^"^, the amplitude of 
fluctuations of the superfluid density becomes of the order of the average 6Ns ~ (Ng); the 
local value of the superfluid density, averaged over the size C,{H), becomes of random sign 
and the system is unstable with respect to spontaneous creation of supercurrents. 

When \H — H^\/H^ ^ G^^, one can neglect the second term in brackets in Eq.ll. 
Rescaling r and A(r) as 

r=^l/Geo,A(r) = ^/(^), (26) 

we obtain a dimensionless equation for f{y) which represents a continuous version of XY 
spin-glass model 



Vlf{y) + j dy'Skiy, y')f{y') = f{y), (27) 

where {6k{y,y')) = and the correlation function 
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{6k{y„ y\)6k{y„ y',)) = 6{y, - y,)6{y\ - y',) ^\ + 5{y, - y2)5{y[ - y'^Vi - v'l)- (2^ 

[yi - y i) 

In this limit, it follows from Eqs. 26-27 that the amplitude of spatial fluctuation of the mod- 
ulus of the order parameter 5A(r) ~ {A{H)) ~ Aq/G is of order of its average. The char- 
acteristic spatial scale of the fluctuations of 5A(r) is of order of Ld. The important feature 
of Eq.27 is that the sign of the second term in Eq.27 fluctuates randomly which corresponds 
to the random sign of the superfluid density. The spontaneously created supercurrents in 
this case have random directions, their typical amplitude is of order of J|^ ~ N^h/G^^o and 
their characteristic scale of spatial correlations is also of order of L^,. The current in Eq.8 
is negligible compared with J'^'^ when l/lso ^ G'^, where Iso ~ VDt^. 

The fact that the sign of A^^ is random is especially clear in the case of a large magnetic 
field, when H — H^ ^ H^G~^. In this case, A(r) can be nonzero only due to existence the 
rare regions, where 6k{x,x') is much larger than the typical value given in Eq.28. Thus, the 
spatial dependence of the modulus of the order parameter has the form of superconducting 
domains embedded in a normal metal. These regions are connected via the Josephson effect. 
We calculate the average critical current of the junctions 



(Jc) oc G — -2 exp --—, 



(29) 



which decays exponentially with the average distance between the superconducting droplets 
Lq. On the other hand, the amplitude of fluctuations of Jc decays only as a power of Lq. 



2 -X 2 



HZ 



iiSJcY) oc ( -^ . (30) 



As a result, the amplitude of the fluctuations in this regime turns out to be larger than the 
average, hence Jc has a random sign. As argued in section 5, such a distribution function 
of exchange interaction is a generic one when the spins are polarized. 

We should emphasis that the Josephson coupling in Eq.30 is derived in the limit when 
Lq is much longer than the coherence length of the superconducting domains. It doesn't 
depend on the value of Aq in the superconducting domains. This is because the effective 
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transmission coefficient of Cooper pairs over a distance of Lq is exponentially small at an 
energy higher than D/L"^ <^ Aq. In contrast, when Lq <^ ,^0) the effective transmission 
coefficient of Cooper pairs is independent of e at an energy smaller than Aq. In this limit, 
Jc oc Aq. 

It is well known [^ that at T = the long range order of the ground state of the 
two-dimensional XY model is destroyed by an arbitrary small concentration of "antifer- 
romagnetic" bounds. As we have mentioned above in the case H <^ H^, regions where 
Ns{r) < 0, exist with small but finite probability. In this case, however, the properties 
of a superconductor are different from the XY model because the supercurrents sponta- 
neously created in these regions are screened by the Meissner effect. Thus at if, T = a 
superconducting ffim should exhibit the conventional long range order. 

This implies that there is a critical magnetic field Hsg < H^ atT = 0, where the system 
undergoes a phase transition from a superconducting state to a superconducting glass state. 
The typical distance of the rare regions is estimated in the next section as L^iSee Eq. 63). At 
Hsg, if should become comparable with the penetration depth, \{H) = X^/dx {H/H^ — HY, 
or 

XiHsc) ~ msG) exp(G^ ^°~f^^ ). (31) 



c 



It yields 

Ao is the zero temperature penetration depth of a bulk superconductor. Eqs. 7,32 show that 
up to a log-factor, the transition between the superconductor and the superconducting glass 
state takes place at the magnetic field when 6 A ~ A(ii). The interval of magnetic fields 
where the system is in the superconducting glass state is indicated in Fig.l. 

The superconducting glass state which arises due to orbital magnetic field effects has 
been considered in numerous papers (See for example [E^-pi]]). The qualitative difference 



between |l29|-pl| and above considered cases is that in the latter case the system exhibits 
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glassy behavior as soon as H > Hd and vortices begin to penetrate the superconductor. Fur- 
thermore, the state we discussed here is also different from Fulde-Ferrell-Larkin-Ovchinnikov 



state [^,Q, which becomes a metal stable state in dirty superconductors when ^p ^ / pi 



IV. OPTIMAL SUPERCONDUCTING DROPLETS: H > H^ 

One of the consequences of the mechanism discussed above is the instability of the spin 
polarized disordered Fermi liquid well above the critical magnetic field. As argued before, 
though the average curvature of the normal metal state {0{r, r') evaluated at A = 0) is 
positive defined, its mesoscopic fluctuations have random signs because of the mesoscopic 
fluctuations of the spin polarization energy. In the regions where the spin polarization energy 
cost to form superconducting pairing state is much lower than the average energy cost, the 
fluctuations of the curvature are of large negative value comparable to its positive average 
such that the normal metal with A = becomes unstable. As a result, above the critical 
field H^, the superconducting pairing correlations are established at mesoscopic scales in the 
different regions in the normal metal and couple with each other via exchange interactions 
of random signs. This argument was present in another early paper by the author [^ . 



In this section, we study the probability to find regions where the superconducting pairing 
states are formed at mesoscopic scales at H > H^. At high magnetic fields in the strong 
spin-orbit scattering limit, the statistics of these pairing states can be studied with the help 
of the generalized Landau-Ginsburg equation, 

A3(r) 



^2/vv .2e \2 H^ H 



A(r) + j iK'iv, r', ff) A(r')<ir' = ^, (33) 



which is valid when H — H^ is small compared with H'^ and when the spatial variation of 



the pairing wave function A(r) over distance ^o is negligible. Here ^o = JD/Aq, A is the 
vector potential of external perpendicular magnetic field. 

Eq.33 is a nonlinear equation in terms of A(r), with a nonlocal (5i^°(r,r') potential 
originating from the oscillations of the wave functions of cooper pairs. Generally speaking. 
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it is qualitatively different from the Scliroedinger equation of an electron in the presence 
of random impurity potentials p6|-|39|. These complications arise naturally in the study of 



the interplay between the mesoscopic effects and the superconductivity and are the generic 
features of strongly correlated mesoscopic systems. In fact, this nonlocal structure of the 
potential in Eq.33 leads to the superconducting glass state. 

At H — H^ ^ H^/G^, the optimal configurations which determine the macroscopic 
properties of the sample turn out to be the superconducting droplets embedded inside the 
disordered Fermi liquid, with the phases of each droplet coupled via random exchange in- 
teraction. Such a configuration can be characterized by three parameters: A), the typical 
size of the droplet, Lj; B). the typical distance between the droplets, L^ ^ Lf] C). the 
typical value of the order parameter inside each droplet. In the following, we will discuss the 
statistics of the mesoscopic pairing states in this regime. In the leading order of (Lj/L^)^ 
the statistical property of the formation of superconducting pairing states at mesoscopic 
scales is similar to that of the impurity band tails p6|-p9|]. 

The calculation of such a probability is closely connected to the evaluation of tails of 
distribution functions of mesoscopic fluctuations |^^-^. However, in the present case. 



(5ii'°(r, r') is determined by the fluctuations integrated over the whole energy spectrum 
instead of single energy level. Thus, we believe it is of a Gaussian form and the statistical 
property of the random potential 6K^{r,r) is determined by its second moment. General 
case is discussed in Appendix C. 

The pairing wave function of the most probable configurations is given as 

M^) = Yl ^aVa{r), / drria{r)r]fs{r) oc Sa/s- (34) 

Note ?7(r) introduced in this way is dimensionless. For such a configuration to have lower 
energy than the normal state, 

/ drdr'A{r)0{r, r') A(r') < 0, (35) 

where 0{r, r') is given by Eq. 19. 
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The total energy of such a configuration consists of cross terms corresponding to the cou- 
phng between different droplets. The coupling between the droplets decays as the distance 
increases. When the size of the droplets is much smaller than the distance between them, 
the typical magnitude of the coupling between different droplets is much smaller than that 
of the coupling within one droplet. We are going to neglect such terms in the estimate of 
the probability of the droplets in the leading order of o{L'l/L'^). 

Thus, to have / droplets in the normal metal, / independent inequalities have to be 
satisfied 



J drUr) ko V + ^^^ Va{r) + J drdr'r]^{r)6K\r, r')r]^{r') 



< 0. (36) 



(We assume there is no perpendicular magnetic field.) Furthermore, we can write down 
the probability to have superconducting pairing states ai H ^ H^ in term of the sum of 
probability to have certain number of droplets 

ViMx)})=J2Pii{r^^}\a = l,...J). (37) 

I 

To simplify the notation, we introduce 

Olg = Co'V^ + ^^^, Km = 6K\v, v'). (38) 

c 

Taking into account Dri{r) = UaDrja, we have 

Pi{{r]^}\a = l,...,l) = Jp{{Km})II^N' Je{-E^ + F^)Dr]^DKM (39) 

where 

Eaiivir)}) = J drr]a{r)OLGVa{r), 

FaiMr)}) = I drdr'r]ir)^KMir, v')r]^ir'), (40) 

and A^ is a normalization constant. P{{Km}) is the distribution function of Km] Df]-, DKm 
represent functional integrals. We use the following equality to transform the step function 
into integrals. 



/U n+oo 

dg^ / dha exp [i/i„ {E^ - Fa - ^a)] • (41) 

Eq. 37 is reduced to -Pi({^a}|ck = 1,...,/) = Tl^Pa- In the Gaussian approximation, the 
statistics of 5K^{y^v') is completely determined by the second moment of the correlation 
function, or 

< Ill^,5K{v,, v[) >= n- iC(r2„ 4; r^.+i, 4+i), (42) 

where C(ri, r'^; r2, r2) is given in Eq.l5. 

In this case, pa can be simplified in a closed form as 

aa= dridr[ / dr2dr2C{ri, r[; r2, r2)ria{ri)ria{r[)r]a{r2)r]a{r2) (43) 



with erfc(a/6) = /^ exp(— a;^/26^)/v27rP^(ia;. One can evaluate the functional integral Drjaijc) 
in the saddle point approximation as long as H — H^ ^ H^/G^. The saddle point equation 
of Eq.43 can be obtained by minimizing the argument of the error function. 

OLG^a(r) + Sy dridr[dr'C{r,r';ri,r[)r]airi)r]c,{r[)r]air') = 

g ^ J drr]a{r)OLGVair) 

J dridr[dr2dr'2C{ri, r[; rs, r'2)r]a{ri)r]a{r[)r]a{r'2)r]a{r2) 

(44) 

The solution of the saddle point equation i]s{r) determines the shape of the optimal droplets. 
To carry out the functional integral of ^a(r), one can expand ri{r) around the saddle point, 

?7(r) = r]s{r) + 6r]{r), 6r]{r) = ^ a„r7„(r) (45) 

n 

where 77„(r) are the eigenstates of the operator r(r, r') generated via second functional 
derivative of the argument in the error function with respect to ri{r) at ri{r) = i]s{r). Our 
final result barely depends on the detailed structure of r(r, r') and we do not give an explicit 
form here. Performing the Gaussian integral of 6ri{r) around the saddle point, taking into 
account the normalization condition, we obtain, 
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E, ^ det'T{r,r') r 

where 

E,= drr]s{r)OLGr]s{r) 

Es/\/2as is the argument of error function in Eq.43 evaluated at ?7(r) = 77s (r). ' indicates the 
exclusion of the zero eigenvalue. The last integral in Eq.46 corresponds to the contribution 
from the zero eigenvalue state, originating from the translation invariance of the saddle point 
equation, with 2 — fold degeneracy 

Lpdoirisir - To) 



i = x,y [^,^. Here Lq is the characteristic length of the droplets determined via the 
normalization condition 



Thus, 



/ [dao] = 72 / dxadya- (50) 

J L/n Jv„ 



The spatial integral is performed only in the region f„ where no other droplets are present. 
Using the following rescaling 



r = yLf, V = Vj^L.\ risir) = r]s 






C(r,r';ri,r;) = ^,%C{y,y'-y,,y[\ Lf=Uj7^f'\ (51) 



we can express E^., (Tg in term of dimensionless T]s{y) 

B = J dyris{y){-Vl + l)risiy),A'^ = J dyidy[dy2dy'2C{yi,y[;y2,y2)Vsiyi)Vs{y'i)Vsiy2)Vs{y2) (52) 
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where B, A? are the dimensionless quantities of order of unity depending on the details of 
Vs{y)- Vs satisfies the dimensionless saddle point equation 



-Vl + l)risiy) + / dyidy[dy'C{y,y';yi,y[)r]s{yi)ris{y[)ri,{y') = 



and at y = oo, ris{y) 
estimate that 



(53) 



0. If i]s{y) is a Gaussian function, B = 2, A = IS/vr^. We also 



cietr(r,r') 
det{0{r,r')) 



l + o( 



G" 



(54) 



Collecting all the results, we have 



Pii{Va}\a = l,--,0 = — ^T2^^f^ 



n I L 



Hr ) 



-I \ I 



A 



(55) 



where V^ /l\ is from the spatial integral in Eq.46, excluding the overlap between different 
droplets. We take into account Lq ~ Lf. It is easy to confirm that the average number 
density of the droplets is 



P 



lElpll 


= ^erfc 


' BG 


(H-H'X^'' 


VEiPi 


Lj 


A 


[ H^ J \ 



(56) 



The distribution function of the amplitude of the order parameter A in a droplet can be 
calculated in a similar way. In this case, the amplitude A is determined by the nonlinear 
term in Eq.33 and the probability to have a superconducting droplet with order parameter 
equal to A is 



A 



P(A) = / -—^P{Km)S{N^{—)^ + E^- F^)dKMd7]a 

J Aq Aq 

where N^ is given as 



(57) 



(58) 



Transforming 6 function into an integral and carrying out the Gaussian integral, we obtain. 



P(A) = N f Dr]air 



2AiV„ 



A 



I ^2na^iMr)}) 



exp[ 



(E^iMr)}) + N^ii- 



A_N2\2 



2a«({r7(r)}) 



(59) 
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The saddle point equation of Eq.59 is similar to Eq.43 except there is an additional nonlinear 
term proportional to A^^. As we will see that the typical A in optimal droplet is much 



smaller than AqJH — H^/H^^, this new term is much smaller than the linear term and can 
be neglected as far as the spatial dependence is concerned. We can use the saddle point 
solution obtained in Eq.43 to evaluate Eq.59, 

, , ,, 2NsA 1 . /.^ . detrir, r') , , 

^ ' Al y/2^ ^^ 2as'det<0{r,r')> ^ ' 



where 



Es + N,^. (61) 



Ns is the corresponding value of Na evaluated at ^^^(r) = ris{r). Substituting the results in 
Eqs.51,52 into Eq.60, we obtain the conditional distribution function of A 



2CGA A' 

-^ exp(-C G -. 

^0 ^0 



P,(A^) = ^-^ eM-C'G'^). (62) 



V. EXCHANGE INTERACTION BETWEEN THE DROPLETS 

The coupling between different droplets deserves special attention. Though the coupling 
between droplets does not affect the probability of finding one droplet, it determines the 
global phase rigidity. The typical distance L^ is of order 

L,erfc-V^[^(^^)V.] (63) 



c 



following Eq.56. It is important that as long as H — He/ He '^ 1, L^ is much less than 
C,c = /exp(— G^), the localization length in the presence of a parallel magnetic field; the 
weak localization effect in this case is small as far as the superconductivity is concerned. 
The typical coupling between a and j3 droplet is determined by 60. Taking into account 
Eq.51, in the limit L^ ^ Lf we obtain the variance of the coupling 

An/L,^2 



'< i^UoAe^Ap J drdr'6K0{r,r%,{r - r„)^,(r' - r^))' > oc ^ (^^) . (64) 



22 



(The coupling depends on Aq in this case because the spectrum in the superconducting 
droplet is gap less.) To get this result, we take into account that the size of the droplet is 
Lf, typical A^ is given by Eq.62 and Ir^ — r^| ~ L^. One the other hand, the average O, 
as shown in Eq.l9 is proportional to 6{r — r'). The average coupling is proportional to the 
overlap integral of the wave functions of two droplets 

/dr,,(r-r.),.(r-r,)«exp(-WL,). (66) 

The variance of the coupling evaluated in Eq.64 is much larger than the average coupling in 
the limit Ld^ Lf. The distribution function of the coupling between different pairing states 
is symmetric with respect to zero and the sign of the coupling between different mesoscopic 
pairing states (droplets) is random. This suggests that the ground state of these coupled 
mesoscopic pairing states will exhibit glassy behavior in this limit. 

The existence of random Josephson coupling in the presence of a parallel magnetic field is 
a consequence of the Pauli spin polarization. This phenomenon exists even without spin orbit 
scattering. Consider for example a granular superconductor, with superconducting grains 
embedded inside a noninteracting disordered metal coupled with each other via Josephson 
coupling. The sign of the Josephson coupling is determined by the total phase of the time 
reversal pairs. In the pure limit, though the sign of the wave function of each electron 
oscillates with a period of the Fermi wave length, the total phase of (p, — p) pair is zero 
because of the exact cancellations of the phases of each electron inside the pair. Therefore 
there is no sign oscillation for Josephson couplings. In the dirty case p is not a good 
quantum number. However the sign of the coupling does not oscillate as a function of 
spatial coordinate because of the time reversal symmetry. As a result, even when the 
distance between the grains is much larger than the mean free path, the sign of the coupling 
is positive defined |^3|]. This is in contrast to RKKY exchange interaction between nuclear 
spins. RKKY coupling exhibits Friedel oscillations with the period of the Fermi wave length 
in the pure case; in the presence of impurity scattering, the phase of Friedel oscillations of 
electron wave functions becomes random. 
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In the presence of a parallel magnetic field, the electrons inside the normal metal become 
polarized. In this case, the electron with spin up has a different kinetic energy as the electron 
with spin down on the Fermi surface because of the Pauli spin polarization. As a result, the 
phase of the electron with the spin up does not cancel with that of the spin down one in the 
presence of Zeeman splitting, and the total phase is equal to Jq Pup ' dr + J^. Pdown ■ dr = 
J(2 drfiBH/vp- Integral is carried out along trajectory C along which electron pairs travel. 
In the pure limit, ~ L^bH/vf with L the distance between two grains. The pairing wave 
function oscillates and develops nodes in its spatial dependence 

2nkT E <- V^^Gf (r, rOG:r^(r, r') = ^^ cos( '^ ~ ^^'^^^ - ^) (66) 

This leads to the sign oscillations of the Josephson coupling with a period Vp/f^'sH, 
which is much longer than the Fermi wave length. The positions of these nodes in the 
spatial dependence of the coupling can be shifted in random directions when impurities are 
present. To estimate these random phase shifts, consider disordered metals with short mean 
free path and L ^ I. The trajectory of electron pairs is a diffusion path with typical length 
L^//. In this case, ~ ^bHL^/D. When L'^/l ^ vf/^bH, is much larger than unity, and 
the sign of the coupling becomes unpredictable for different impurity configurations. In this 
limit, the Josephson coupling averaged over impurity configurations is exponentially small 



exp{—\/2L/JD/fiBH) while the typical amplitude of the coupling decays as L''^. Therefore 
when the magnetic field increases, only the position of the maximum of the distribution 
function moves towards zero while the width of the distribution function barely changes. 
This results in the superconducting glass state. Note that in principle the charging effect 
inside the grain will also lead to the superconducting glass phase as suggested in a recent 
experiment ||1|. However in the metallic limit when the tunneling conductance between the 
grain and the normal metal is much larger than e'^/h, charging effect should be negligible 
and only the mesoscopic mechanism discussed in this paper is relevant. 

In this section, we find that BCS order parameter is determined by mesoscopic fluctu- 
ations of physical quantities. Short range mesoscopic fluctuations are responsible for the 
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presence of optimal superconducting droplets while long wave length fluctuations lead to 
frustrations. Close or above the mean fleld critical points, the inhomogeneous supercon- 
ducting states are described by a nonlocal Landau -Ginsburg theory. 

VI. MESOSCOPIC SAMPLE 

For a flnite system of size L there are, in principle, many critical flelds H^. Linearizing 
Eq.33 with respect to A, neglecting the gradient term and using perturbation theory we 
have 

6H, = H,- i/o = Hl^ J 6K%r, r', H'^ + 6H,)drdr' (67) 

To derive Eq.67 we have taken into account: 1. The relative amplitude of fluctuations of the 
critical field is smaller than its average {{6Hc)'^) ^ {H^Y. 2. The sample size is smaller than 
the coherence length L < ^{H) and A(r) is spatially uniform. Eq.67 reflects the fact that 
the magnetic fleld acts on the system in two ways: a) It suppresses superconductivity, b) 
It changes the mesoscopic fluctuations of parameters of the normal metal and the quantity 
J 6K^{r,r',H)dr'dr is a random function of 6Hc. Therefore, generally speaking, at a given 
T, Eq.67 can have an inflnite number of solutions, which means that the iJ-dependence 
of the critical temperature Tc{H) exhibits reentrant superconductor- metal transitions as a 
function of H. Qualitatively the picture of the reentrant transitions is very similar to that 
which takes place in the case of magnetic fleld induced orbital effects [|14| . To characterize 



the random quantity 6Hc, we study the statistics of s, which is the right hand side of Eq.67. 
Straightforward calculation of its variance following Eq.l5 yields 

(To =< s^ >= / dridr2C{ri, r[; r2, r'a) (68) 

Its distribution function in the Gaussian limit reads as 

Following Eq.67, the distribution of H^ is 
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P\H,) = J5{H, - Hi - Hls{H,))P{s{H,))ds{H,) = _exp(- 



1 1 , ,H?,-H,^^ 1 



2aoH, ^^ ^ H^ ' 2^0 '■ 

(70) 



In deriving the second line we use that 



<^.>^0. (71) 



It is obvious following Eq.70 that the variance of H^. is 



(72) 



Eq.72 gives the interval of the magnetic field near H^ where the reentrance takes place with 
a probability of order of unity. The probability for a sample in a superconducting state at 
H can be estimated as 

V,{H) = r°^ P%H)dH = ^^ erfc{G ^^~^^ ). (73) 

When spin-orbit scatterings are weak, AqTsq ^ 1, the conventional theory leads to the 
conclusion that the superconductor-normal metal transition is of first order with the critical 
magnetic field Hcc 0,0- In this case the spin polarization in the superconducting phase 
is zero. The average spin polarization energy of a normal metal sample of size L and its 
mesoscopic fluctuations are of order dL'^i'o^fiBHY and dLC,oh'Q{fiBH)^G~^ , respectively. As 
a result, a finite superconducting sample exhibits first order normal metal -superconductor 
reentrant transitions in the interval of magnetic field of order 6H = Hcc^q/GL in the vicinity 
of the critical field. 

In the case of two dimensional superconducting film, the fluctuations of both polarization 
energy of the normal metal and the condensation energy of the superconducting phase should 
lead to a nonuniform state, similar to the case Aqt^o ^ 1. The theory of this phenomenon 
at AqTso ^ 1 is, however, more difficult. In this case a domains of normal phase within a 
bulk superconductor (or a superconducting domain in normal metal) has the surface energy 
of order of dLj:,C,oi'oAQ, where L^, is the domain size. This energy is much larger than the 
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typical energy associated with mesoscopic fluctuations in Eq.36. Thus the probability of the 
occurrence of such domains is small even at the critical point. 

It is worth emphasising that qualitatively, the case AqTso ^ 1 is not different from the 
case AqTso ^ 1 for in both cases the superconducting glass solutions survive at T = and 
H > H^. Especially, for a quasi ID thin stripe with the width W <^ L^, the surface energy 
becomes independent of the size of the domain while the mesoscopic fluctuations of spin 
polarization energy is proportional to \fT^. This situation is similar to the strong spin 
orbit scattering limit discussed before. 

VII. CONCLUSION 

We show the existence of a novel superconducting glass phase in disordered thin films 
in Clogston limit. The statistics of mesoscopic pairing states in the superconducting glass 
phase is universal and determined only by the sheet conductance. It is a direct consequency 
of Wigner-Dyson statistics of single particle energy spectrum. 

This allows us to distinguish the mechanism discussed in this paper and the effect of 
inhomogeneity of impurity concentration, or classical pinning effect on vortex lattices dis- 
cussed in [Q. First of all, in the present case, the magnetic field couples only with spins 
and the wave functions are real (as far as the impurity averaged condensate wave function 
is concerned); the time reversal symmetry is broken spontaneously. For classical pinning 
effects on vortex lattices, the time reversal symmetry is broken by the applied perpendicular 
magnetic field. More over, fluctuations of local quantities like mean free path can lead to 
inhomogeneous states but do not lead to spontaneous time reversal symmetry breaking. The 
glass state discussed in this paper is due to random signs of long range exchange interaction, 
which is purely of mesoscopic nature. Finally, the response of the state discussed here is 
determined universally by Thouless energy of the size of the coherence length and the re- 
sponse of a pinned vortex glass depends very much on the range and strength of the classical 
pinning potential. For amorphous films where the impurity potential is perfectly screened 
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and in the absence of granularities, the classical pinning effect is weak; the mesoscopic effects 



dominate in this limit. Amorphous thin films like Zn |^, Mo-Ge |^6[, Pb ||5^, In-InO |^ 



Bi [^] have been subjects of extensive studies. 

Though the transport properties of such a superconducting glass state are poorly un- 
derstood, it shares all the features a glass state has: hysteresis, stretched relaxation time. 
Another experimental consequence of random sign of Ns{r) we like to mention is, follow- 



ing to Ref. [^], at H > Hsg and at a finite temperature the system exhibits the negative 
magnetoresistance with respect to the component of the magnetic field perpendicular to the 
film. 

For a finite sample, when gate voltages are applied, the mesoscopic fluctuations in Eq.67 
start to oscillate. This causes reentrant superconductor-metal phase transitions, similar to 
the magnetic field induced reentrant transitions. Such transitions should manifest them- 
selves in the gate voltage finger print experiment: the conductance as a function of gate 
voltages exhibits sample specific fluctuations, with amplitude equal to the normal sample 
conductance. The conductance fluctuation can much exceed the value of UCF due to the 



attractive interaction! Q . 

The other possibility to study the mesoscopic superconductor is to bring the supercon- 
ducting state adiabatically along a closed trajectory in a parameter space via applying gate 
voltages(for thin films like Bi, the chemical potential can be varied by 20 percent.). Adiabatic 
charge transport across a boundary of the system per period in the presence of periodically 
changing external perturbations is connected with a geometric phase, as first pointed out 
by D. Thouless 0]. Recently, this idea was applied to normal metal mesoscopic systems 
where quantum chaos is fully developed; the charge transport is determined by the amount 
of "flux" of a topological field which threads the area enclosed by a closed trajectory in 
the parameter space pT] , p2| . In a normal metal mesoscopic sample, such a topological field 



was shown to be determined by the sensitivity of the quantum chaos to external pertur- 
bations, which is a random quantity. In the case of superconductors, the geometric phase 
will be determined by the compressibility of the superfluid density because excess electronic 
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density created by gate voltages can be carried away only by coherent motions of the con- 
densate. The mesoscopic fluctuations of superfluid density are "more compressible" than 
the electronic density itself, i.e. 

/i is the chemical potential. The existence of mesoscopic fluctuations of condensate wave 
function or superfluid density manifests itself in a geometric phase. This problem will be 
addressed elsewhere. 

The question whether or not the quantum fluctuations of the phase of the order parameter 
destroy the superconducting glass state at T = and large H is still open. We would like 
to mention here, that a similar question was addressed in many papers in the context of the 



disorder driven superconductor-insulator transition |[53|-|55| and metallic spin glasses with 
dissipation [^ . 



Strictly speaking, at arbitrary finite temperatures T > 0, a superconducting film doesn't 
possess a superconducting phase rigidity. In two dimensional case, due to screening, the 
interaction energy between vortices decays as a power law rather than logarithmically. This 
leads to a finite concentration of unbounded vortices with the correlation function of the 



phase of the order parameter decaying exponentially at large distances 1^7 . 

However, in real experiment situations, the London penetration length can be compa- 
rable or longer than the sample size. Furthermore, the exchange interaction decays as a 
power law function, 1/r^, as shown in Eqs. 15,28. The typical energy of a domain of size 
L diverges logarithmically as L — > oo, suggesting that there could be a finite temperature 
phase transition between the superconducting glass phase and a normal phase. In Fig.l, we 
plot a dashed line which separates these two phases, the existence of which needs further 
investigation. On the other hand the two dimensional XY model with short range random 
exchange interaction is known not to exhibit a phase transition between the paramagnetic 
and the spin- glass phases ||5B |. 



The frustration which leads to the novel superconducting glass phase is due to the ex- 
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istence of nodes in the condensate wave function when spins are polarized, as emphasized 
in section 5. For a d-wave superconductor, nodes exist even in the absence of an external 
magnetic field; naturally one can ask whether a d-wave superconductor can be free of frus- 
tration when disordered. We are not aware of work on this subject and believe the answer 
to this question is also critical to the understanding of the density of states at the Fermi 
surface in a disordered d-wave superconductor. 

The author acknowledges useful discussions with B. Altshuler, C Biagini, D. Huse, S. 
Kivelson, I. Smolyarenko, B. Spivak, N. Wingreen. This work is supported by ARO under 
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VIII. APPENDIX 

A. Fluctuations of spin polarization energy 

Diagrams in Fig. 2d yield 



liSTP \'i\ — /'0-TrT"^2Y^r /3i/3i^/3iQ!2 ^-Q2-a2r^"2/32 _/32/32/^/32Q!i ^-ai-ai^Qi/31 
\\'J^p) I — K'^T^J- ) Z^Fz <-^ei+e2'^2 '-^ei+e2'^2 ^ei+e2^z '-"ei+e2 

irrPiPi r)l3ia2 _-Q2-a2 n"2/32 -ftfe n^^ai _-ai-ai 7-)Q!i/3i 1 /yp-N 

+ '^2 ^ei+e2"z ^<i\+<i'i"z ^ei+e2"z ^e-i+e2\- \'^) 

0""° = a, a = ±1 is for spin index; e = (2n + l)iTkT is Matsubara frequency. J2 represents 
summations over Matsubara frequency, momentum and spin index. 
Following Dyson equation in Fig. 2c, 



2 VI + ^q' + Tsoif^sHy \iu\ + Dq2 + r,-i^ 

D°^(q2) = i[ ^ + 1 1 (76) 

^^^' 2V| + Dq2^ |cc;| + Z}q2 + r,-iJ ^ ^ 

The sensitivity to the change of pair potential 6 A is given by Eq.75 with C, D replaced via 
6C,6D 

^r+.2 (q')c^2T (o)c^2T WDTJ% (q')^2^? (o)c^2T (o)^J.t2 (q') 

C7+.2 (q')^2".? (o)c^2T (o)c7;re2 (q')cie? (o)c^2l'^« (o)c:,t2 (q') , (77) 

following Fig.2e. 

B. Fluctuations of superfluid density 

The correlation function of the mesoscopic fluctuations of the superfluid density consists 
of two terms. First term is given in Fig.2g 
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eie2 

(78) 
while the second part of the contribution in Fig.2h 

{6u>6jj>6iri - r2)(5(r; - r'2) + (5,,v(5,,/(ri - r^)(5(r2 - r[)}N^{H) x 

C'2r^ (0)^2^4"^ (0)^2^^ (0)^2!^^ (0)Dj;7t. (q')I^Zf?^^ 

(79) 

C. Pa when 50 is non- Gaussian 

In general, statistics of 60 is characterized by the following correlators, 

< n^50(r„ r[) >= C2™(r„ r^l^ = 1, .., 2m), (80) 

Eq. 43 then is transformed into Pi{{ria}\(^ = 1, •••,0 = ^aPa, where 

/O r+00 

dga / dhaDr]a{r) exp{iha{La - go)) 
-00 J —00 

1.2m « 

^^ J llk=i,...,2mVa{rk)v»{r'k)drkdr'f.C2m{{rk,r'k}\k = l,...,2m) 

= li"^5l7r~r^I^^^/(^"'^^) / ^k=i,...,2mVa{rk)Va{r'i,)drkdr'^.C2m{{rk,r'k}\k= l,...,2m). (81) 

The main contribution is from the saddle point where 
6p{Mr)}) 



6ri{r) 



OLGV{r) + V{{r]{r)})=0 



1 9^™ f 

li";iZl7r~r^^^^/(^"'^^) / ^k=i,...,2mVa{rk)Vair'k)drkdr'^5{r - rk)C2mi{rk,r'k}\k = l,...,2m). 

(82) 
In the Gaussian limit, Eq.82 yields Eq.43. 
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FIGURES 
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FIG. 1. Suggested phase diagram in Clogston limit when AqTsq <^ 1; Insert is the 
phase diagram of Chandrasekhar-Clogston theory. A superconducting glass phase appears at 
HsG = H^ — 5Hc with 5Hc ~ H^/G'^lnG. The thick line along T = axis represents the 
superconducting glass phase discussed in the paper; the dashed line stands for a possible finite 
temperature phase transition between the superconducting glass phase and a normal phase. 
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FIG. 2. a). Dyson equation for diffusons and Cooperons. Solid lines correspond to electron 
Green functions in a metal; Dashed lines are for impurity scatterings preserving the spin while 
the wavy ones for spin orbit scatterings, b) .Vertex correction of pair potentials, c). Hikami 
Box. d),e). Diagrams for the fluctuations of spin polarization energy, f) Diagrams representing 
the correlation function {6K^{ri,r2)SK^{r3,r4)). g),h). Diagrams representing the correlation 
function of supercurrent densities (Js(r)Js(r')). Solid wavy lines in g). represent the correlation 
function {5A{ri)6A{r2)) given in Eq. 17. 
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